Mesoscopic spin Hall effect along a potential step in graphene 
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We consider a straight one-dimensional potential step created across a graphene flake. Charge and 
spin transport through such a potential step are studied in the presence of both intrinsic and extrinsic 
(Rashba) spin-orbit coupling (SOC). At normal incidence electrons are completely reflected when 
the Rashba interaction (with strength A^) is dominant whereas they are perfectly transmitted if the 
two types of SOC are exactly balanced. At normal incidence, the transmission probability of the step 
is thus controlled continuously from to 1 by tuning the ratio of the two types of SOC. Besides the 
transport of charge in the direction normal to the barrier, we show the existence of a spin transport 
along the barrier. The magnitude of the spin Hall current is determined by a subtle interplay 
between the height of the potential step and the position of Fermi energy. It is demonstrated that 
contributions from inter-band matrix elements and evanescent modes are dominant in spin transport. 
Moreover, in the case of vanishing extrinsic SOC (Xr = 0), each channel carries a conserved spin 
current, in contrast to the general case of a finite Xr, in which only integrated spin current is a 
conserved quantity. Finally, we provide a quasi-classical picture of the charge and spin transport by 
imaging flow lines over the entire sample and Veselago lensing (negative refraction) in the case of a 
p — n junction. 



I. INTRODUCTION 

Within the last decade the role of SOC in pure crystals 
has been fundamentaly reconsidered thereby providing a 
host of novel effects [1-4] and phases [5-12] with potential 
applications in spintronics. 

First it was predicted that SOC may generate a trans- 
verse spin current in response to an applied electric field 
both in hole and electron doped semiconductors [1, 2]. 
This so-called intrinsic spin Hall effect was promptly re- 
ported by the detection of the related spin accumulation 
at the boundaries of GaAs samples [3, 4]. 

More recently SOC was shown to lead to a topological 
phase of electronic matter when combined with a partic- 
ular band inversion property [5-14]. In two-dimensional 
systems this so-called Quantum Spin Hall (QSH) state is 
characterized by metallic edge states surrounding an in- 
sulating bulk. Those edge states show up in the absence 
of magnetic field and realize a time- reversal invariant ver- 
sion of the chiral edge states of the Integer Quantum 
Hall states. Initially the QSH state was introduced in 
the Kane-Mele model of graphene which describes non- 
interacting electrons on honeycomb lattice with both in- 
trinsic and extrinsic SOC [5, 6]. In the Kane-Mele model, 
intrinsic SOC conserves the z-componcnt the real spin 
(s z ) and generates a topological mass term. In contrast 
the extrinsic Rashba-type SOC breaks s z conservation 
by mixing s z = +1/2 (spin f) and s z = —1/2 (spin J,) 
spin components. When the intrinsic SOC dominates 
the Rashba one, the bulk excitations are gapped and a 
pair of gapless states counter-propagate along an edge 
of the sample. Nevertheless the weakness of the SOC 
induced gap [15-17] makes the realization of the QSH 



state extremely difficult in graphene. In contrast sim- 
ilar QSH edge states have been predicted [7] and soon 
after reported in transport experiments, performed in 
HgTe/CdTe quantum wells [8, 13, 14]. 

In a recent paper [18], we have shown that charge 
transport through a potential step allows to investigate 
the relative strength of intrinsic and extrinsic SOCs. 
In the absence of SOC, the transmission of electro- 
static potential steps have been extensively studied in 
graphene both experimentaly [19-24] and theoretically 
[25-27] with the purpose of providing a condensed-matter 
implementation of the relativistic Klein tunneling. In 
practise such potential steps can be induced either by a 
distant gate [19-24] or by metallic contacts [28, 29]. 

In this paper we detail supplementary aspects of the 
charge transport properties outlined in [18]. In addition 
we emphasize here the spin transport near an electro- 
static potential step. As a main result, an interfacial 
spin Hall effect is predicted, which we call hereafter meso- 
scopic spin Hall effect (MSHE). In the MSHE the spin 
current flows in the direction transverse to the applied 
electric field, and is mainly localized at the vicinity of 
the step. Owing to this nonuniform spin current distri- 
bution, the MSHE differs from the spin Hall effect in ho- 
mogeneous doped semiconductors [1-4] and in graphene 
[30]. 

The paper is organized as follows. In Sec. II, the 
charge and spin current operators for the Kane-Mele 
model are derived in presence of intrinsic and Rashba SO 
coupling. It is shown in section III that quantum averages 
of these operators consist of direct and cross (interfer- 
ence) terms when computed within a generic scattering 
state. We discuss charge transport in Section IV. As a 
main prediction of this paper, the spin transport along 
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the interface (transverse to the applied electric field) is 
described thoroughly in the section V while possible ex- 
perimental detection is also discussed therein. Finally, 
we provide a quasi-classical picture of the charge and 
spin transport by imaging flow lines on the entire sample 
and Veselago lensing (or negative refraction) at the p — n 
junction. 



where indices a, (3 = ±1 specify the band. The dispersion 
relation of the band a(3 reads 



E a p(p) = a^p 2 + (A + f3\ R ) 2 + p\ R , 



(3) 



and the corresponding wavefunction ipafip( x ) = i x \ a P) p 
can be expressed as 



II. FORMALISM: CHARGE AND SPIN 
CURRENTS 

We consider a graphene monolayer in the presence of 
both intrinsic and extrinsic SOC . The corresponding 
charge and spin current operators are constructed in the 
framework of the Kane-Mele model of graphene [5, 6]. 
The quantum averages of the charge and spin currents are 
derived both for propagative and for evanescent quasipar- 
ticlcs. 



A. Kane-Mele model 

The low-energy Kane-Mele model [5, 6] is defined by 
the Hamiltonian Hkm = Hq + Hso + Hr with 

Hso = -^Aa z r z s z ip, 
H R = ^ \r{<j v s x - a x T z s y )ip. (1) 

Here cri,Ti,Si (i — x,y,z) are the Pauli matrices associ- 
ated with the lattice isospin (A and B sites of the hon- 
eycomb lattice), the valley- isospin (K and K' points of 
the reciprocal space), and the real electronic spin, respec- 
tively. 

The kinetic Hamiltonian Hq describes graphene in the 
absence of any spin-orbit interaction. The intrinsic spin- 
orbit effect (with coupling constant A) is described by 
Hso which preserves all the symmetries of the problem 
and further conserves the component s z of the electronic 
spin. In contrast the Rashba contribution Hr (with cou- 
pling constant Xr) explicitely breaks the conservation of 
Sz ■ 

The full Kane-Mele Hamiltonian Hkm acts on 8- 
spinors. Nevertheless in this paper we will only consider 
intravalley scattering caused by electrostatic potential 
steps. Hence we shall focus on the Zf-valley in the fol- 
lowing analysis and simply substitute (t z = 1) in Eq.(l). 

The resulting Hamiltonian H^m consists of a 4 x 4 ma- 
trix acting on a spinor of the form [ipAf, V'sfi i^b\\i 
where f, J, stands for real spin. 

In the homogeneous case (in the absence of spatially 
varying potential), the momentum p = (p x ,p y ) is a good 

quantum number. The single-valley Hamiltonian H^ M 
is diagonalized by the eigenstates \a/3) p as 



Px ~ iPy 

E af} (p) + A 
-ip{E a p{p) + A) 
-ip{p x +ip y ) 



(4) 



with 



P 



WL. 



(5) 



These wavefunctions arc normalized to represent a unit 
probability within a rectangular graphene flake of length 
L (along the x-direction) and width W (along the y- 
direction). 

In the presence of a potential step or a sample edge, 
evanescent states appears. Their spectrum and wave 
functions arc given, respectively, by Eqs. (3,4,5) with 
typically an imaginary p x . Note also that p 2 = p 2 — 

\p x \ 2 ■ Boundedness of the wave function does not al- 
low evanescent modes to exist in an ideally infinite sys- 
tem. They are, nevertheless, ubiquitous in any hetero- 
junctions (near potential steps) or even at a sample ter- 
mination. The helical edge states, on the other hand, oc- 
cur only at the boundary between two topologically dis- 
tinguishable phases each characterized by opposing Zi 
indices. In the case of Z 2 topological insulator, Zi in- 
dex distinguishes trivial (even number of Kramer's pairs) 
vs. non-trivial (odd number of Kramer pairs) insulating 
states. 



B. Charge and spin current operators 

We construct here the charge and spin current opera- 
tors associated with the single- valley Kane-Mele Hamil- 
tonian H^^. We simply write the continuity equations 
for the charge and the s z -componcnt of the electronic 
spin. The charge and spin density operators (p c (x) and 
p Sz (x) respectively) can be written as, 



p c {x) = -ei'^x^ix), 



(0) 



using the four component field operators i){x) and ip'(x). 
The charge current operator J c is determined such that 
it satisfies the continuity equation for charge: 



H 



(K) 
KM 



\a/3) p = E a p(p)\aP) p , 



(2) 



dp c 
dt 



divJ c = 0. 



(8) 
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Using the equation of motion 

dip 

i— = [—itr ■ V + \ R {o y s x - <j x Sy) - Aa z s z )} -0, (9) 

it is readily verified that the operator defined as, 

J c = —etp^crip, 

satisfies Eq. (8). This expression can be alternatively 
obtained from the usual definition of the charge current, 



dH 



{K) 

KM 



op 

and by noting that the only momentum-dependent part 
of H^ M is Ho = i^p ■ a-ip, (v F = 1)- 

We repeat the same procedure by introducing the spin 
current operator for s z -spin component as 



J Sz = -^(TS z tp. 



(10) 



However, in the presence of Rashba SO coupling, the sec- 
ond term on the r.h.s. of Eq. (9) does not commute with 
s z , leading to a source term in the continuity equation 
for spin: 



dp Sz 
dt 



div J, 



)j>. (11) 



The source term Xr(u x s x + o~ y s y ) describes the spin 
torque due to Rashba spin orbit interaction. Violation of 
the continuity equation leads in many cases to the ambi- 
guity in the definition of spin current density, reflecting 
correctly the fact that spin is not conserved. However, 
as far as spin transport in our junction problem is con- 
cerned, we will see later in Sec. V that the spin current 
density defined as Eq. (fO) yields a conserved spin cur- 
rent when integrated over the incident angle <f>. Such a 
spin current satisfies the continuity equation globally, i.e. 
in the sense of 



PF 



2ir 



dt 



+ div(J s J =0, 



(12) 



where (...) denotes averaging over the scattering states. 
The spin current density (10) becomes locally conserved 
in the absence of Rashba SO coupling: \r = 0. 



C. Charge and spin current carried by an 
eigenstate \a/3) p 

We now evaluate the quantum average of the charge 
and Sz-spin currents in the eigenstate \aj3) p defined by 
Eq. (4). When both p x and p y are real, this state cor- 
responds to a propagating plane wave and carries the 
charge current density 



(Jc) 



afip 



-AeA 2 a JE aP (p) + A) p 



(13) 



which is collinear to the momentum p. This current con- 
sists in equal contributions from the channels s z = 1 and 
s z = — 1. Moreover these contributions correspond to 
opposite spin currents leading to a cancellation of spin 
current for the s z -spin component. This statement is 
confirmed by direct calculation of the quantum average 
of the s z -spin current J s , , 



(JsJ 



0, 



(14) 



in the bulk eigenstate \a/3) p . 

In presence of a potential step or at a sample edge, 
evanescent states, described by an imaginary p x , become 
possible. If we assume a semi-infinite graphene plane 
extending over the half-plane (x > 0), the charge current 
carried by such an evanescent wave reads 



{Jc)ct/3 P 



0. 



teAl p {E af) (p) + A) pyt 



(15) 



The current is localized near the interface and flows along 
the y— axis. Moreover the net charge transport vanishes 
when the sum over p y is performed. 

Finally the average of the spin current in an evanescent 
state, 



( J ^)al3p _ > 

{ J i) a p P = - 2A l^ E ^iP)+^) 1^ 



(17) 

-2|p.|* ( 18 ) 



shares the characteristics of the charge current except 
that spin current will not vanish upon p y integration (see 
section V). 



III. SCATTERING STATES FOR A POTENTIAL 

STEP 

We consider a potential step in a graphene monolayer 
as shown in Fig. 1 and construct the corresponding scat- 
tering states. We evaluate the quantum averages of the 
charge and spin current densities in a given scattering 
state. Owing to the spin-orbit coupling (i.e. multiband 
character of the Kane-Mclc model) , the structure of those 
averages is more complicated than the average currents 
in a pure eigenstate |a/3) p . Indeed a single incident elec- 
tron generates two transmitted electronic waves. There- 
fore the averaged currents contain direct terms involv- 
ing only one kind of quasiparticles, and crossed terms 
describing coherent interferences between the two trans- 
mitted waves. 



A. The junction model 

We consider an electrostatic gate creating a potential 
step V(x) which is slowly varying on the scale of the 
atomic lattice. Moreover we assume an ideally pure sys- 
tem which can be approached in experiments with sus- 
pended devices. These conditions ensure that intcrvallcy 
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FIG. 1: Potential step V(x). 



scattering can be neglected and that the junction is de- 
scribed by the Hamiltonian 



H = H$ [ +tfV(x)a s il>, 



(19) 



where <7q arL d sq represent the identity in lattice isospin 
and real spin space, respectively. Besides, we further as- 
sume that the step is sharp on the scale of the Fermi 
wavelength. In this sense we may represent the potential 
V(x) by an abrupt step defined by 



V(x) 



(x < 0) 
V (x>0) 



(20) 



Note that we also assume a straight interface with trans- 
lational invariance along the y direction (no roughness 
along the interface x = 0). 



B. Scattering states 

Typical scattering state with energy E and transverse 
momentum p y 



*E, Py (x) 



ipi(x) + rip R (x) + r ev ip ev (x) (x < 0) 
t + ^> + (x) +t-ip-(x) = iPt{x) (x > 0) ' 

(21) 

can be constructed in terms of the bulk eigenstates Eq. 
(4). It is assumed here that the incident electron is a pure 
| + — ) p state. Hence the incident and reflected waves read 



i>i, R {x) = A(E,p) 



±Px ~ iPy 

E + A 
i{E + A) 
i(±p x + ip y ) 



i(±p x x+p y y) 



where ^(resp. — p x ) corresponds to the incident -0/ 
(resp. reflected 4>r) wave, and 



1/A(p) 



p 2 + {e + ay 



WL. 



The x-component of the momentum p x = p x (E,p y ) is 
the positive solution of E^ {PxiPy) = E. There is also 



an evanescent state defined in the half-plane x < 0, 



1pev{x) = A{E,p) 



-IPy 

A 



—IK 

E - 
i(E + A) 
+ ipy) 



D KX+ip y y 



and localized near the interface x = 0. The value of k = 

n(E,p y ) is set by the positive solution of E^ (iK,p y ) = 

E 

Within the half-plane x > 0, the scattering state con- 
sists of two transmitted waves with opposite symmetry 
f3 = + and ft = — , described by the spinors 



^(sc) = A{E,p) 



Pfix ~ IPy 

E + A 

-i0 (E + A) 
-i/3(p/3x + ip y ) 



j(Pf3:cX+p y y) 



where the x-components of the momentum pp x = 
pp x (E,p y ) are obtained by solving E a p{pfj x , p y ) = E-V . 
Note that pp x can be either real or imaginary. When pp x 
is real, 4>p stands for a propagating mode and the sign 
of p x is chosen such that the group velocity is positive, 
i.e. correctly describes an outgoing wave. When pp x 
is imaginary, ipp is an evanescent mode and the sign of 
the imaginary part of pp x is fixed by the requirement of 
wavefunction boundncss at x — )■ oo. 

It is sometimes more convenient to specify ^E,p (x) by 
E and an incident angle <fi, instead of E and the trans- 
verse momentum p y : 



P.i 



\p\ COS ( 



Py 



p| sine 



(22) 



Finally the four scattering amplitudes r, r ev , t+ and t_ 
are uniquely determined by solving the continuity condi- 
tion at x = [18]: 

^E, Py (x = 0- , y) = V E , Py (x = 0+ , y) (23) 

for given E, Vq and p y (or (j>). 

C. Direct and crossed terms 

We have seen that the scattering state ^ E,p y (x) takes 
the form of a superposition of two bulk states with oppo- 
site band symmetry /3. The expectation value of current 
density in such scattering state 

J c = -eip^o-ipT 

= - e (|< + | 2 VV0+ + M 2 <AW- 



(24) 



has naturally two types of contributions: direct and 
crossed terms. The direct terms (proportional to \t + \ 2 
and \t~\ 2 ) arc similar to the ones addressed in the pre- 
vious section, see Eqs. (13,14,15,16,17,18). In particular 
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such direct terms were shown to carry no net s 2 -spin cur- 
rent when associated with a propagative wave (Eq.(14)) 
whereas evanescent waves carry a finite spin current 
(Eq.(17)). In contrast the crossed terms 2$l[t* i _t-^ + (T'i/j-} 
always contribute to spin transport regardless of the na- 
ture of the two interfering transmitted waves. 

In the following we focus on the charge and spin trans- 
port associated with these latter crossed terms which 
mix the transmitted waves ip+ an d ip- altogether. The 
crossed charge current is proportional to the expression 



= 2A + A-(E + A) 



P-x - P +x 



J{P-x-P+Jx 



(25) 



A. Pseudo- reflection symmetry 

The continuum limit of the Kane-Mele model, defined 
as in Eq. (1) has a pseudo-reflection symmetry (PRS) 
operating within each valley. PRS with respect to the 
x-axis: y — > —y is expressed as U — o x s y V, where V 
represents a parity operator in two spatial dimensions, 
i.e.,p — > p' = VpV~ x = (p x , —p y )- The eigenstate of the 
Kane-Mele Hamiltonian, Eq. (1), is also an eigenstate of 
PRS operator U, since [H, U] = 0, with an eigenvalue, 
—f3, i.e., Ui/j a /3 P (x) = —f3tjj a pp>(x'). This can be checked 
explicitly using Eq. (4). We demonstrate that PRS is a 
convenient tool for clarifying the dependence of reflection 
and transmission coefficients on incident angles. 



where E = E — Vb , while p +x and p- x were defined in the 
previous subsection. These crossed terms yield a charge 
current along the y-axis. In contrast to the direct terms, 
the crossed current has a spacial dependence upon the 
coordinate x (dependence upon y is forbidden by trans- 
lational invariance). The total current is therefore diver- 
genceless. 

We now consider the spin current. The crossed terms 
are proportional to: 



= 2A + A-(E + A) 



P-x +P +X 

2 P y 



(26) 



The spatial distribution of the spin current is oscillatory 
if both transmitted states (ip+ and ip-) & re propagative. 
If one of the transmitted states is propagative while the 
other is evanescent, the spin current distribution shows 
damped oscillations. 

The x-component of Eq. (26) is generally finite and 
has a x-dependence. The y-component of Eq. (26) is 
also finite, but does not contribution to the divergence of 
spin current density. Therefore, 



divJ, 



dx 



(27) 



remains finite due to the cross term. Recall that a spin 
current-density is generally not a conserved quantity, see 
Eq. (11). In contrast the contribution from direct terms 
to the spin current is divergenceless. A similar circum- 
stance also occurs in a more conventional semiconductor- 
based spin Hall system [31]. 



IV. CHARGE TRANSPORT 

We consider the charge transport across an electro- 
static potential step, in the presence of spin-orbit cou- 
pling. The charge conductance is directly determined by 
the transmission probability whose energy dependences 
are investigated thoroughly in this section. 



B. Transmission at normal incidence 

The continuity equation at x — reads, 

ipr + ri>R + r ev ^ ev = + t-V>-! (28) 

where y-componcnt of momentum is +p y , determining 
the incident angle 4> for a given p x . Applying PRS oper- 
ator U from the left to the above equation, one finds, 



ipj + rip' R - r ev ip' ev = -t + ip' + + t-tp'_ 



(29) 



This can be regarded as the continuity equation for 
momentum— p yi where ip' = ^p\ Py= - Py . Comparing the 
two equations, Eqs. (28) and (29), one can convince one- 
self that reflection and transmission coefficient are either 
an even or an odd function of the incident angle: 



fey j 



t' = -t + , t'_ = t- 



(30) 



At normal incidence: p y = 0, two states with an opposing 
PRS eigenvalue, — /3 are orthogonal to each other at the 
spinor level, i.e., 



^a/J( P:c ,0)( )</W'(^0)(0) OC <V 



(31) 



for arbitrary p x and p' x . Consequently, the continuity 
condition Eq. (23) is decoupled to two equations: 



ipl + ripn = t-ip-, 



(32) 
(33) 



as far as p y = at the interface x = 0. The reflec- 
tion coefficient r is determined only by the transmitted 
state with j3 — — , therefore the normal incident charge 
transport is independent of states of the symmetry differ- 
ent from the incident state. Namely, the transition from 
P = — to P = + is impossible due to different symmetry. 

Solving Eqs. (32) and (33) yields the reflection coeffi- 
cient 



X-Y 
X + Y' 



(34) 
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0.0 u 1 1 1 1 1 1 1 

0.5 1.0 1.5 2.0 2.5 3.0 



FIG. 2: Transmission probability T(p y = 0) at the normal 
incidence as a function of Vb/A and E/A at Xr/A = 1.2 
(upper), and the dependence of T(p y = 0) on Xr/A for a 
fixed value of (E/A, Vb/A) = (0.5,2) (lower). In the upper 
panel, each region corresponds to a finite range of T(p y = 
0) indicated on the panel. In the lower panel, each curve 
corresponds to different values of Vb/A. 



where X = pp(E — Vq + A) with pp being the Fermi mo- 
mentum in the incident side, Y = p- x (E+A). Reflection 
probability is given by R = 1 — |r| 2 [18]. Imaginary p_ x 

which corresponds to the evanescent state | ) leads 

to perfect reflection R = 1 even if the other transmitted 
state | — h) is propagating. This perfect reflection oc- 
curs provided A — 2Xr < E — Vq < —A in the dominant 
Rashba case (Xr > A). By contrast, the perfect trans- 
mission r = occurs at the phase boundary A_r = A. 
Therefore the crossover from perfect reflection to perfect 
transmission occurs with tuning Rashba SO coupling. 

The upper panel of Fig. 2 shows the normal incident 
transmission probability T as a function of Vq and E. 
The (black) diagonal strip is the region of perfect reflec- 
tion, which occurs due to dominant Rashba effect. With 
fixing Vq,E and decreasing Xr/A, the reflection proba- 
bility decreases as shown in the lower panel of Fig. 2. 
In the figure all the lines go through the point being 
Xr/A = l,T|0 = o = 1, namely perfect transmission al- 



ways occurs, because the Dirac cone appears again at 
the balanced Rashba and intrinsic SO case (Xr/A = 1). 
Further decreasing Xr/A, the transmission probability 
decreases and reaches to zero at Xr/A = (1 — E/A)/2 = 
0.25. Crossover from perfect reflection to perfect trans- 
mission occurs due to the competition between Rashba 
and intrinsic SO in the normal incident case. 



C. Charge conductance 

Let us consider the rectangular geometry (Fig. 1) with 
infinite aspect ratio W/L, W and L being respectively the 
width and length. We have considered so far the charge 
current carried by a single scattering channel with def- 
inite p y . In the absence of disorder, the channels are 
independent and the charge conductance G c (in units of 
e 2 /h) is simply the sum of the single- channel transmis- 
sions T(p y ) over all possible transverse momenta p y , or 
equivalent ly over all incidence angles: 

C ^r®T W . (35, 

We evaluate Eq. (35) explicitly by substituting T(p y ) = 
1 — \r(p y )\ 2 , where the reflection amplitude r(p v ) follows 
from the continuity condition (23). Note T ^ |M 2 + |*-| 2 
because of potential difference Vq between the incident 
and transmitted side. The obtained charge conductance 
is shown in Fig. 3 as a function of Vq for different values 
of Xr/A. The curves exhibit specific features depending 
on the value of Xr/A. We leave further inspection of 
such behaviors to Sec. Ill C. 

Fig. 3 shows several conductance curves as a function 
of Vb/A for different values of Xr. The incident energy is 
set to be E/A = 0.5. The conductance is normalized by 
Go = Wpp/i: where pf is the Fermi wavevector in the 
incident side. 



1. Semimetalic phase (Xr > A) 

When Xr = 1.7 (Fig. 3(a), red curve), the system is in 
the semimctallic phase. The charge neutrality (particle- 
hole symmetric) point on the transmitted side is located 
at Vb/A = E/A + 1 = 1.5. The conductance vanishes 
at this point. Above this value of Vo, the conductance 
shows a singularity at Vb/A = E/A - 1 + 2X R /A = 2.9 
(discontinuity in the first derivative). Here, the Fermi 

energy touches the lowest energy band: | ) on the 

transmitted side. Above this value, the number of energy 
bands contributing to the conductance is doubled. The 
conductance shows an abrupt increase of slope (a kink) 
at this point (see also Appendix). 

Below the neutrality point (Vb/A = 1.5), the con- 
ductance shows a peak at Vq = 0, then turns to a 
slow and monotonic decrease as Vb is decreased. This 
feature does not seem to resemble its behavior above 
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FIG. 3: Upper panel: charge conductance G c of the potential 
step (height Vb) normalized by Go = Wpf/tt (E/A = 0.5). 
Lower panel: Fano factor F. Red, blue and green curves 
represent a typical behavior of charge conductance, respec- 
tively, in the semi-metallic phase (Xr/A = 1.7), in the topo- 
logical gap phase (Xr/ A = 0.6) and at the phase boundary 
(An/A = 1). 



the neutrality point. Decreasing Vq from Vq = 0, the 
Fermi energy touches the highest energy band: | + +) at 
Vq/A = E/A - 1 - 2X R /A = -3.9. But the conductance 
curve does not show any singularity here. 

This asymmetric behavior is a fingerprint of the unique 
band structure of Kane-Mele model. At Vq/A = 2.9, 

the | ) band touching the Fermi energy has the same 

band index /3 = — 1, consequently the same symmetry 
as the incident state: | H — ). Therefore, as soon as this 
state becomes available for transport, transmission oc- 
curs, typically, at the normal incidence, leading to an 
abrupt increase of the conductance. On the other hand, 
at Vb/A = —3.9, the | + +} band touching the Fermi sur- 
face has the symmetry opposite to that of the incident 
state. Therefore, no transmission occurs at the normal 
incidence via the | + +} band, and the conductance curve 
bears only a gradual change, (see also Appendix A). 

Note also that perfect reflection in the semimetallic 
phase is limited to normal incidence, and the conductance 
takes generally a finite value due to contribution from 
non-zero incident angle, (f> =/= 0, to the r.h.s. of Eq. (35). 
In the case of bilaycr graphene, the conductance curve 
shows similar features, with two shoulders only on the 
V Q > side. 



2. Topological gap phase (Xr < A) 

When A^ = 0.6 (Fig. 3(c), blue curve), the system 
bears a band gap. Naturally, the conductance vanishes 
identically in the gap region: —A < E — Vq < A — 

2X R , i.e., 0.7 <V /A< 1.5. The | } band is always 

activated for transport in the pn-rcgimc: Vq/A > E / A + 
1 = 1.5, showing no singular behavior in the conductance 
curve. 



3. At the phase boundary (Xr = A) 

At the phase boundary (Fig. 3(b), green curve), two 
bands with /3 = — arc combined, and a pair of lin- 
early dispersing energy bands, i.e., a Dirac cone ap- 
pears. The conductance curve is similar to that of 
monolayer graphene without spin-orbit interaction. [27] 
The conductance vanishes at the charge-neutrality point 
Vq = 1.5, and shows again no singularity on the pn-side 



D. Fano factor 

Fano factor F associated with the potential step is ob- 
tained [27] from the transmission probability T as 



F 



pf 



d Py T{l - T) 



PF 



PF 



dp y T. 



(36) 



PF 



The lower panel of Fig. 3 shows the Fano factor as a func- 
tion of Vq for different values of Xr. One of the specific 
features is that the Fano factor shows, independently of 
the value of A^, a peak structure at Vq = E + A = 1.5, 
with a maximal value of F (= 1). This is, of course, 
partly related to the fact that the conductance vanishes 
at this point. However, recall also that in monolayer 
graphene without SO interaction, the Fano factor re- 
mains structureless at the corresponding point Vq = E 
[27], with a value of F « 0.1. Such characteristic sup- 
pression of shot noise is spoiled by the SO effects, even 
when Xr = A, where Dirac cone reappears. A zero of 
the conductance leads, trivially, to a peak in the Fano 
factor. Very contrastingly, the charge conductance is not 
qualitatively affected by the SO effects. 

The Rashba dominant regime exhibits some anomalous 
features: a cusp in the conductance at Vq/A = 2.9 for 
Xr/A = 1.7 (case (d) in Fig. 3) and an enhancement of 
the Fano factor, in the regime of perfect reflection. This 
corresponds to 1.5 < V /A < 2.9 for X R /A = 1.7. In the 
region of perfect reflection, evanescent modes are dom- 
inant in transport at a finite incident angle. The Fano 
factor is enhanced by such evanescent modes. The Fano 
factor increases with the increase of Rashba SO coupling 
in this regime of perfect reflection. In the balanced case, 
Xr = A (Fig. 3(e)), the conductance curves do not differ 
significantly from the case of no SO effects (Xr = A = 0). 
In contrast, when intrinsic SO coupling dominates the 
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0.4 0.9 1.4 1.9 2.4 2.9 

FIG. 4: Charge conductance (upper) vs. Fano factor (lower) 
as a function of Rashba SO coupling. Incident energy is set 
to be E/A = 0.5. Each curve corresponds to different values 
of Vb/A. 



Rashba term, i.e., in Fig 3(f), the conductance vanishes 
within a finite range of Vq, corresponding to the gap. The 
Fano factor is not well-defined, therefore, not plotted in 
this regime of Vq . 

In the large and small enough potential region 
|Vo/A| 3> 1 and for arbitrary Xr, the Fano factor takes 
a value F « 0.1, which is roughly the same as the Fano 
factor without SO effects, since the transmitted state is 
free from evanescent modes in this region. 



E. Crossover effects in conductance and Fano 
factor 

The charge conductance and the Fano factor are shown 
in Fig. 4 as a function of Rashba SO coupling Xr/A. 
At normal incidence, we have seen a crossover from per- 
fect reflection to perfect transmission in the pn-regime: 
Vq/A > 1.5 (curves (i) and (ii) in the lower panel of Fig. 
2). As for transport properties, perfect transmission at 
normal incidence is replaced by a broad maximum near 
A_r = A (see cases (g) and (h) in the upper panel of 
Fig. 4). The maximum always appears in the pn- regime. 
However, unlike the normal incident case, the maximum 
of charge conductance is not unity and is not precisely 
located at Xr = A as a consequence of angular integra- 
tion. 



The charge conductance becomes smaller with the de- 
crease of the Vq on the pn-side (e.g., cases (g) and (h) 
in the upper panel of Fig. 4), and actually vanishes at 
Vq/A = 1.5 (not plotted). This is because the number of 
propagating states on the transmitted side reduces with 
the decrease of the Vq (and vanish at Vb/A = 1.5). A 
cusp of conductance appears also appears in this regime 

at the band edge of | ). When Xr/A is larger than 

this value, the Fermi energy intersects with only (3 = + 
band. As a result, perfect reflection occurs at the normal 
incidence, leading to smaller values of conductance. 

On the other side, i.e., in the nn-regime: Vb/A < 1.5 
(e.g., cases (i) and (j)), maximum does not appear, be- 
cause perfect reflection docs not occur at the normal in- 
cidence (see (hi) and (iv) in the lower panel of Fig. 2). 
When Vb/A = 1 (case (i)) , charge conductance vanishes 
in Xr/A < (l-E/A+V /A)/2 = 0.75, where the system 
is gapped on the transmitted side. 

Dependence of Fano factor on Rashba SO coupling is 
shown in the lower part of Fig. 4. Interestingly, the lower 
panel looks almost the upside down image of the upper 
panel. 

V. SPIN TRANSPORT 

Here we investigate the spin transport generated by an 
electrostatic potential step in the presence of spin-orbit 
effects. The potential step splits the graphene sample in 
two pieces characterized by distinct carrier densities. The 
spin is drifted along the interface and the spin current is 
therefore transverse to the applied electric field. This 
spin Hall effect (SHE) is a mesoscopic analog of the bulk 
spin Hall effect occurring in homogeneous electron or hole 
doped semiconductors [1, 2]. The present effect requires 
both spin-orbit coupling and a step in the carrier density. 
The spin Hall current localized in the vicinity of interface 
appears also in 2DEG with Rashba spin-orbit interaction 
[32]. 

A. Local spin current density 

The spin flows along the y-axis, namely along the in- 
terface x = defined by the potential step. The total 
spin current density, 

= *T IZ ff*^K^(*), (37) 

results from the summation of the single-channel currents 
J| (a;) over all possible transverse modes labeled by their 
momenta p y . This local current is maximal near the in- 
terface and decays when the distance x is increased (Fig. 
5). Note that the spacial dependence of spin current den- 
sity j'Jf (x) is determined by a subtle interplay between 
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V /A = 5 
V /A = 2 
V /A= 1.4 




FIG. 5: Spatial distribution of spin current density. Xr/A = 
2 (0.5) in the upper (lower) panel. Spin Hall current is local- 
ized in the vicinity of the junction, jo = A/(2irL). 



the direct and cross terms. When both transmitted waves 
are propagative, the (singlechannel) crossed spin current 
oscillates function of the distance x from the in- 

terface. Of course, these oscillations arc damped when 
integrated over all possible transverse momenta p y (see 
Fig. 5) but this decay is a lonf ranged power law rather 
than an exponential decay. 

When one or two transmitted wave(s) is/are evanes- 
cent, even the singlechannel current decays exponentially 
when the distance x from the interface is increased. As a 
result, the total spin current density decays very abruptly 
with x (see Fig. 5). 



B. Spin conservation 

The same situation occurs for the x-componcnt of spin 
current which comes only from the cross term, i.e. from 
Eq. (26). Moreover this contribution, due to spin torque 
in the presence of Rashba SO coupling, vanishes after in- 
tegration over incident angle. The spin transport occurs, 
therefore, only in the direction parallel to the interface. 
Besides, the continuity of spin current-density is recov- 
ered after this angular averaging: 



dp y 
~2T 



dt 



div (J Sz ) 



1.5 
1.0 
0.5 
0.0 
-0.5 
-1.0 
-1.5 



-X s /A = 2 
-A«/A = 1 
— \ R /A = 0.5 















y /A 




The 



; spin current density, defined as Eq. (10), is thus 
considered to be a conserved quantity in the 



sense of 



1.60 



FIG. 6: Spin transport in the Kane-Mele p-n junction for dif- 
ferent values of Ah. The p-n junction induces s z spin current 
in the direction parallel to the interface. Such a spin Hall 
current is plotted at x = as a function of Vq. 



Eq. (12), and free from the usual issue of defining a spin 
current. 



C. Results 

Fig. 5 shows the spatial distribution of spin current 
density parallel to the interface, i.e., y-component of spin 
current is plotted as a function of coordinate x normal to 
the interface for different values of Rashba SO coupling: 
Xr/A = 2 (dominant Rashba phase) for the upper panel 
of Fig. 5, and A/j/A = 0.5 (topological gap phase) for 
the lower panel Fig. 5. 

In the upper panel, 1.5 < Vo/A < 3.5 corresponds to 
the regime of perfect reflection, and the absolute value 
of spin current is large compared with other cases. The 
spin current is localized in the vicinity of interface, explic- 
itly manifesting that spin is carried by evanescent modes 
(localized in the ^-direction but propagating in the y- 
direction). It also shows a damped oscillatory behavior 
for a larger value of Vo such as Vo/A = 5 in the upper 
panel. Such damped oscillation is an incarnation of cross 
terms between evanescent and propagating modes. The 
lower panel of Fig. 5 corresponds to the topological gap 
phase, in which the spin current takes a large negative 
value when Fermi energy is in the gap on the transmitted 
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side. 

Fig. 6 shows spin Hall current at x = 0. Let us com- 
pare it with the charge conductance shown in Fig. 3, 
both represented as a function of Vo- The two curves 
show indeed quite contrasting behaviors. At Vo = (in 
the absence of a junction) the charge conductance show a 
maximum (peak). In contrast, the spin current vanishes 
at Vq = 0, reflecting the fact that here, spin transport is 
a mesoscopic effect due to the presence of interface. 

Both magnitude and sign of the spin current is tuned 
by Vq . The direction of spin current is opposite between 
inter-band tunneling (Vo > E + A) and intra-band tun- 
neling (Vq < E + A) cases (compare the two panels of 
Fig. 6). [39] The latter includes also the case of metal- 
insulator junction. We can confirm this explicitly from 
Eqs. (18) and (26) with E E = E - V . The direct 
term of spin Hall current is proportional to E — Vq + A, 
therefore the direction changes near Vq = E + A. The 
sign of the crossed term is opposite from the direct terms 
and enhanced in the vicinity of Dirac point Vq = E + A, 
resulting finite spin current at Vq = E + A. 

In the dominant Rashba phase (Xr > A), two 
quadratic bands touch at the neutrality point E = —A. 
On the transmitted side, this corresponds to Vq = E + 
A = V n . At this value of Vq, the charge conductance 
vanishes. The spin conductance changes its sign in the 
vicinity of neutrality point, but remains finite precisely 
on that point (see the lower panel of Fig. 6). This is again 
due to the evanescent modes. The spin current takes a 
large positive value above the neutrality point: Vq > V n . 
In the topological gap phase (Xr < A), the charge con- 
ductances vanishes in the gap: E—A < Vo < E+A = V n , 
whereas the spin current is enhanced in the gap, taking 
a large negative value. 

The upper panel of Fig. 7 reveals two different natures 
of mesoscopic spin Hall effect by studying, separately, 
the A = and = cases. Fig. 7 shows that the 
mesoscopic spin Hall current flows actually in the absence 
of topological mass term: A = 0. The spin Hall current 
is enhanced in the regime of perfect reflection: E < Vq < 
E + 2\ R (0.5 < Vq/Xr < 2.5 in Fig. 7). This corresponds 
to the regime of Vq above the neutrality point, where the 
spin conductance takes a large positive value. The lower 
panel of Fig. 7 shows, on contrary, J| as a function of Vo 
in the absence of Rashba SOC: Xr ■ — 0. Clearly, the spin 
current is enhanced, when Fermi energy is in the gap on 
the transmitted side: E - A < V Q < E + A = V n (0.5 < 
Vq/A < 2.5 in Fig. 7), i.e., in the situation of metal- 
insulator junction. In the region of Vq < and Vq > 2E, 
the spin current vanishes because only propagating mode 
appears. The spin degeneracy remains in the absence of 
Rashba SOC, as a result, the crossed term also vanishes. 

On the transmitted side, the enhancement of spin Hall 
current thus occurs for two reasons: (i) perfect reflection 
in the dominant Rashba phase (on the Vo > V n side), 
(ii) due to the topological gap (on the Vq < V n side). In 
the two cases, evanescent modes play a dominant role 
in the solution of scattering problem at the junction. 



1.0 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 

V /A 



FIG. 7: Mesoscopic spin Hall effect in the absence of topolog- 
ical mass term: A = (Rashba SO coupling: Xr = 0) in the 
upper (lower) panel. Spin Hall currents at x — are plotted as 
a function of Vo for E/X R = 0.5 with A = 0, jo = X R /(2ixL) 
and E/A = 1.5 with \ R = 0, jo = A/(2ttL). 

The enhancement of spin Hall current along the interface 
thus has two different natures, both related to evanescent 
modes. Depending on which side of the neutrality point 
Fermi energy on the transmitted side is, the enhanced 
spin current flows in the opposite directions. 



VI. ELECTRON VESELAGO LENSING 

The p — n junction in graphene is expected to serve as 
an electronic version of " Veselago lens" [33] . The charge 
current distribution: j c (x) = (j*(x),jv(x)), can be used 
for imaging the electronic flow around the potential step 
[34]. 

Let us imagine an electronic wave packet emitted from 
a point x = (—a, 0) with a > 0. If the wave packet has 
a center-of-mass momentum p = (p x ,Py) with p y /p x = 
tan (f>, then the wave packet will be incident at the p — 
n junction (located at x = 0) at y = atan</>. Let us 
consider the trajectory of this wave packet after it goes 
through the pra-junction. 

Here, instead of following explicitly the dynamics of 
such a wave packet, we calculate directly the stationary 
charge current distribution on the transmitted side, using 
Eq. (24). Then, we consider stream lines of the vector 
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FIG. 8: Refraction of the electron beams by the potential step in the topological gap phase (A_r/A = 0.5). Left panel (a): 
intra-band tunneling (Vo/A = —2). Right panel (b): inter-band tunneling (Vo/A = 2). Spatial coordinates: x = xA/(hvF), 
y = t/A/(fii)j?). Fermi energy: E/A = 0.5. 





FIG. 9: Electron Veselago lens — semimetallic phase (A_r/A = 2). Spatial coordinates: x = xA/(Jivf), y = yA/(hvp). Fermi 
energy: E/A = 0.5. The height of potential step Vb is chosen such that (a) left panel: Vo/A — 2, (b) right panel: Vo/A = 3.5, 
both corresponding to the inter-band tunneling case. Case (b) corresponds to the opening of the lowest energy channel. 



field j c (x). Once j c (x) is known, the following differen- 
tial equation: 



dx 



(38) 



determines the locus of a stream line under a given 
boundary condition, say, y(x — 0) = atan</). Equa- 
tion (38) determines, in turn, the trajectory of the wave 
packet emitted from x = (— a,0). As we have seen in 
Eqs. (13), (14), (15), (16), (25), j% has no y-dependence, 
whereas jjf has no spatial dependence due to charge con- 
servation. Therefore, with a given boundary condition, 
Eq. (38) can be trivially integrated to give, 



y(x) = atancj) + 



1 



dx'jy(x'). 



(39) 



Repeating the same procedure for different values of in- 
cident angle 4>, one can draw a set of stream lines visual- 
izing the vector field j c (x). Focusing of such stream lines 
can be regarded as an electronic version of optical lens. 
Figs. 8 and 9 demonstrate such electron lens realized at 



the Kane-Mele p — n junction. A color code specifies the 
strength of current flow along each trajectory. 



A. Topological gap phase 

In the topogical gap phase, the refraction proper- 
ties are quite similar to the ones studied previously for 
graphene p — n junctions in the absence of spin-orbit in- 
teraction [33]. Indeed the y-component of the current 
density changes its sign on crossing the p — n junction, 
thereby realizing the negative or Veselago-like electronic 
refraction. In presence of SOC, the system does not 
longer show perfect focusing (Fig. 8.b). 

In contrast for a n — n junction (intraband transmis- 
sion) the y-component of the current has the same sign 
on both sides of the junction, indicating that the refrac- 
tive index is positive. As a result, the outgoing electron 
beam is divergent (Fig. 8. a). 
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FIG. 10: Detailed plots of the shade region of Fig. 9 (b): 
Xr/A — 2, Vb/A = 3.5. Other parameters are also the same. 



B. Semimetallic phase 

More interesting are the refraction properties of the 
semimetallic phase. Indeed the evanescent modes mani- 
fest themselves by the bending of the electronic rays on 
the transmitted side (Fig. 9). Moreover no transmission 
is allowed at the normal incidence which yields a shade 
area behind the origin O(0, 0) (Fig. 9). 

Figure 9 shows Vcsclago-like electron lens in the 
semimetallic phase (Xr/A = 2) for two distinct steps V$: 
(a) V /A = 2 (left panel), (b) Vb/A = 3.5 (right panel), 
both corresponding to a p — n junction (inter-band tun- 
neling) . 

In case (b), the Fermi energy is touching the top of 

)-band, which happens when Vq = E — A + 2Xr. 

As the Fermi energy approaches the the top of | )- 

band, the cross term between propagating (due to | — 

+)-band) and evanescent (due to | )-band) modes 

plays a significant role. The meandering stream lines in 
case (b) is a consequence of such " cross term transport" . 
It should be also noticed that the shade area is much 
pronounced in case (b). Indeed, one can observe that 
the refractive index turns to be positive in the vicinity of 
the shade area. Due to perfect reflection at the normal 
incidence, the x-component of current density is virtually 
zero for a small incident angle (j>. As a result, the electron 
beam is strongly refracted, with divergent j y /j x on the 
transmitted side (as <f> — > 0). This leads to the formation 
of shade area (see Fig. 10 for detailed plots). 

Study of electron lens behavior thus reveals rich 
mesoscopic transport properties of the Kane-Mele in- 
junction. It should be underlined that this unique meso- 
scopic transport is carried out by the evanescent modes 
and the cross terms. The former enables transport along 
the edge even when the Fermi level is in the gap on the 
transmitted side [35]. 

Finally the strinking difference between the behaviors 
of the topological gap phase (Fig. 8) and the semimetallic 
phase (9) might be observed by scanning probe measure- 
ments similar to those successfully implemented on top 



of ballistic two-dimensional electron gases [36, 37]. 



VII. CONCLUSION 

We have studied theoretically charge and spin trans- 
port at a potential step (both n — n and p — n junc- 
tions) within the Kane-Mele of graphene. We have high- 
lighted the role of reflection symmetry associated with 
the band index (3 in the crossover from perfect reflec- 
tion to transmission while tuning the Rashba coupling 
Xr. We have also computed experimentally measurable 
quantities such as conductance and Fano factor. 

Due to the multiband character of the model, one inci- 
dent electron yields two distinct transmitted quasiparti- 
cles. The spin Hall current, which is mainly localized in 
the vicinity of interface, results from the superposition of 
two types of contributions: (i) direct terms involving one 
kind of transmitted quasiparticlcs and (ii) cross terms de- 
scribing interferences between the two kinds of transmit- 
ted quasiparticlcs (Sec V). The direct terms were shown 
to carry no net s z -spin current when associated with a 
propagative wave whereas evanescent waves carry a fi- 
nite spin current. In contrast the crossed terms always 
contribute to spin transport regardless of the nature of 
the waves. The interplay between those direct and cross 
terms is also important for charge transport (Sec. IV). 

Moreover the electronic flow exhibits a large variety 
of patterns (Sec. VI). In particular a dominant Rashba 
SOC (semimetallic phase) leads to curved rays owing 
to the presence of evanescent states whereas rays are 
still straight for dominant intrinsic SOC (topological gap 
phase). Note that in a monolayer graphene without SOC, 
stream lines are straight and refracted only at the inter- 
face. In principle, it should be possible to identify those 
contrasted shapes by scanning a charged tip above the 
graphene flake as it was performed for two-dimensional 
electron gases in GaAs heterostructures [36, 37]. Fi- 
nally detecting a clear fingerprint of the role of SOC in 
transport measurements in graphene seems not impos- 
sible but difficult, since the magnitude of SOC is small 
in graphene, at most, on the order of ~ IK. [38]. An 
alternative way to probe such unique transport charac- 
teristics of p — n junction may be to use materials with 
stronger SOC, such as HgTe/CdTe heterostructures. 
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FIG. 11: Contribution of | ) band to charge conductance. 

The charge conductance G- (G+) calculated from j_ is 
shown in the upper (lower) panel. System's parameters are 
chosen such that X R /A - 1.7, E/A = 0.5. 



Appendix: Asymptotic Behavior of Conductance 

The curve of Xr/A = 1.7 in Fig. 3 illustrates 
the characteristic features of charge conductance in the 
semimetallic phase: X R > 1. The conductance curve 
shows, say, at A^ = 1.7, a kink structure at Vq = 

E — A + 2Xr = 2.9, on opening of the | )-channel 

to transmission. The purpose of this appendix is to es- 
timate the asymptotic behavior of G c in the vicinity of 
this singularity. 

In Sec. IV C wc estimated G c , by substituting T((f>) = 
1 — \r(cf>)\ 2 with r(cf>) determined by the continuity condi- 
tion (23), into the Landauer formula (35). Here, to reveal 
the nature of singularity at Vq = E — A + 2A/?, we ex- 
tract the contribution of | )-band to G c , and analyze 

its asymptotic behavior in the vicinity of the singular- 
ity (see Fig. 11). Let us introduce SV, the height of 
potential barrier measured from the singularity: 



When SV < 0, the | )-band gives no contribution 

to the charge current on the transmitted side, since the 
state is evanescent. When SV > 0, a propagating mode 
becomes possible with the momentum 



hvpq. 



V2(A ft - A)SV - (hv F ksh 



(A.2) 



obtained from the energy conservation: E a f}(pp x ,p y ) = 
E — Vq, provided \4>\ < 4> m , where the critical angle (j> m 
is defined by, 

(^ F fcsin0 m ) 2 = 2{X R - A)5V. (A.3) 

In the limit of 6V — > 0, this reduces to, 

<f> m ~ y/2(X R - A)SV/(hv F k). (A.4) 

The charge current transmitted to the lowest energy band 
is given by 



evF 2\t- \ fvupq- 
WL E -V + \r 



(A.5) 



from Eq. (13). Integrating over all incident angles, one 
finds 



86V 



20 m j_(e F ,O) 



E + A 



E + X R V E- A + 2A, 



(A.6) 



SV = Vq — E + A — 2Xi 



(A.l) 



The charge current transmitted to the | )-band thus 

shows a linear uprise (proportional to SV) when SV > 

0. As a result, the charge conductance shows an abrupt 
increase of slope at Vb = E + A. Note that t- (E, 0)^=0 
thanks to the same symmetry, i.e., the same /3 (= — 1), 
of the | )-band as that of the incident energy band: 

! + ->• 

Similarly, wc investigated the asymptotic behavior of 
G c in the vicinity of the opening of | + +)-channel at 
V = E — A — 2X R . However, since the | -I- +)-band (fi- 
nal state) has the opposite symmetry (opposite (3) to the 
initial state: | H — ). Due to this mismatch of symmetry, 
transmission coefficient t + ((f> = 0) vanishes even when 
Vq < E — A — 2 A_r. As a result, the leading order contri- 
bution to the charge conductance starts at second order, 

1. e., oc SV 2 (see Fig. 11, right panel). 
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